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Abstract 


One of the primary requirements for X-33 is that it be capable of flying autonomously. That is, on- 
board computers must be capable of commanding the entire flight from launch to landing, including cases 
where a single engine failure abort occurs. Guidance algorithms meeting these requirements have been 
tested in simulation and have been coded into prototype flight software. These algorithms must be 
sufficiently robust to account for vehicle and environmental dispersions, and must issue commands that 
result in the vehicle operating within all constraints. Continual tests of these algorithms (and 
modifications as necessary) will occur over the next year as the X-33 nears its first flight. This paper 
describes the algorithms in use for X-33 ascent, transition, and entry flight, as well as for the powered phase 
of PowerPack-out (PPO) aborts (equivalent in thrust impact to losing an engine). All following discussion 
refers to these phases of flight when discussing guidance. The paper includes some trajectory results and 
results of dispersion analysis 


Introduction 

In June of 1996, Lockheed Martin SkunkWorks was 
selected for the Phase II contract of the design and 
development of the X-33 vehicle. Marshall Space 
Flight Center's Flight Mechanics, Guidance, 
Navigation and Control Systems Branch was selected by 
LMSW to develop the ascent, transition, entry and abort 
guidance algorithms for the X-33. These algorithms 
have been developed and documented and are undergoing 
simulation testing. Pending large changes in the 
vehicle design or subsystem behavior, the guidance 
algorithms architecture is complete. 

The development of the guidance algorithms was 
challenging as they were being designed, coded, and 
tested in concurrence with the vehicle design and landing 


site selection phases. The ambitious schedule of the X- 
33 project demanded the co-development of on-board 
software during the hardware design phase. Guidance 
algorithms were required to handle all of the conceived 
mission profiles and single-engine failure abort 
missions. The X-33 is a sub-orbital, single-stage, 
autonomous, re-usable launch vehicle. It’s purpose is 
to demonstrate the technologies necessary to develop the 
Reusable Launch Vehicle (RLV), eventually replacing 
the Space Shuttle. 

The X-33 has the capability to launch vertically, 
achieve Mach 11, shutdown it's LOX-LH2 powered 
linear aerospike engines, and coast to a horizontal 
landing. The X-33 is scheduled to launch from Edwards 
Air Force Base (EAFB) in July, 1999. The primary 
landing site is at Michael Army Air Field (MAAF) in 
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Utah, approximately 360 nmi downrangc of the launch 
site. 

There arc a series of test flights to MAAF, each with 
specific goals to test the aerothermal protection system. 
To properly guide the vehicle through all of the mission 
profiles, the guidance and control system was developed 
as generically as possible. Through the use of pre- 
mission I-Ioads, the vehicle has the capability to fly the 
numerous mission profiles and single-engine failure 
aborts corresponding to each of these nominal missions. 
This paper will provide an overview’ of the ascent, 
transition, and entry guidance algorithms for the X-33. 
Also to be discussed is the use of the Performance 
Monitor software for the single-engine abort scenarios. 

Ascent Guidance 

The guidance and control software resides on the Flight 
Manager (FM) processor on the vehicle. Ascent 
guidance is executed at 1 Hz. Fig. 1 is a top-level flow 
diagram of the ascent guidance algorithm. 

The Ascent_Guide function accepts the navigation state 
data and computes the expected relative velocity at the 
end of the guidance cycle using the current inertial 
velocity and acceleration vectors. The openloop 
function is executed from launch to a pre-determined 
time in the trajectory (tgr). The openloop function 
performs a table lookup of the steering attitude 
commands, throttle, and mixture ratio as a function of 
the computed relative velocity. In this manner, the 
guidance is always leading the control system. 

The openloop function is followed by the acceleration 
tracking function. This function modifies the throttle 
command from openloop to track a nominal acceleration 
profile. This function reduces the Main Engine Cut-off 
(MECO) dispersions from uncertainties in engine 
performance and vehicle aerodynamics. 

If the mission warrants the use of closed-loop ascent 
guidance, the function, IGM (Iterative Guidance Mode) 
is executed. IGM is the linear tangent steering law is 
the same as that used for the Saturn V (Iterative 
Guidance Mode-IGM, Reference 2). IGM has been 
slightly modified to accept relative targeting data 
necessary for the sub-orbital X-33 flights. The sub- 
function, igmguid, accepts targeting data computed in 
IGM and performs the attitude command calculations. 



Figure 1 

X-33 Ascent Guidance Flow Diagram 

Unfortunately, the X-33 flights to MAAF do not have 
sufficient time in near-vacuum flight to w’arrant the use 
of IGM. Whether or not a closed- loop vacuum scheme 
is actually used on X-33, it will be needed for the 
follow-on VentureStar. 1 IGM, while available for use, 
is not currently in use on any of the X-33 flight 
simulations. The primary reason for this is that the 
Performance Monitor provides a closed-loop procedure 
that is better for X-33 flight than IGM. This will be 
discussed later in this paper. 

If an engine has failed during flight, the ascent guidance 
is responsible for computing the throttle level necessary 
to avoid cavitating the LOX pump. The LOX pump 
must maintain a minimum axial acceleration level 
(NPSP - Net Positive Suction Pressure) to prevent 
cavitation. This required acceleration level decreases as 
the throttle level decreases. While the engine would 
like a low throttle level, guidance requires a high 


1 Other guidance and control procedures needed for 
VentureStar, such as for automated rendezvous and dock, 
are also being worked at MSEC, but are not part of the 
scope of the X-33 effort. 
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throttle level as one engine (thus one half of the 
expected acceleration) has failed. The function, 
NPSP_lhrottIe, computes the highest possible throttle 
level which satisfies the acceleration requirements of the 
engine. 

Throughout the ascent portion of flight, the current 
attitude is computed in the Ascent_Guide function given 
the attitude quaternion from navigation. The ascent 
guidance functions employ a launch inertial plumbline 
coordinate system. The attitude command is computed 
internally in the plumbline system and converted to a 
unit attitude quaternion and sent to the control system. 
Given the current and commanded attitude, the ascent 
guidance function also computes the body attitude rates 
necessary to achieve the commanded attitude by the next 
execution of Ascent_Guide. 


representing desired angle of attack and velocity bank 
angle magnitude. The appropriate bank sign to use is 
computed at transition initiation and is a simple 
function of the MECO velocity heading and the desired 
velocity heading (refer to Notation section) 


s gn(cr. w/ ) = -sgn(i// ~V d ) (1) 

See Figure 2 for definitions of heading angles. The 
bank sign thus computed is used throughout transition 
and subsequently used to initialize entry guidance bank 
command sign. In order to provide smooth bank and 
angle of attack commands, the output <7 cmd of the filter 

represented by the following equation is passed as bank 
command to the attitude control system 


The throttle command is also sent to the ascent control 
system where it is broken down into four chamber 
pressure commands. Ascent flight control has engine 
differential throttling and aerosurface deflection 
capability to control the attitude of the vehicle. 

The mixture ratio command from ascent guidance is 
sent to the Propellant Utilization function. This 
function monitors the propellant levels in each tank and 
modifies the mixture ratio command to reduce the 
propellant residuals at MECO. 

The X-33 has unique engine shutdown requirements. 
Traditionally, the engine shutdown command is issued 
from ascent guidance when the target MECO condition 
is achieved. A flight performance reserve in the form of 
propellant is 

loaded on-board to guarantee that the target will be met 
despite vehicle and environmental uncertainties. 
However, the X-33 has grown in weight to the point 
that the landing gear will not support much residual 
propellant. Thus, the engine begins it's own shutdown 
procedure based on LOX depletion sensors located in the 
LOX feed lines. Ascent guidance continues to be 
executed until the engine shuts itself down. This effects 
MECO state dispersions in that propellant loading 
uncertainties extend or reduce the burntime of the 
engines. 

Transition Guidance 


A™, + 2 $(0 n <7 cmd + or n (a cmd - a ti ) = 0 ( 2 ) 

and analogously for angle of attack 

a«r + 2go)„cc cmd + - a ti ) = 0 (3) 

The filter states (X cmd and G cmd are initialized during 

the first pass through transition guidance by setting the 
commanded values equal to the current estimated bank 
and angle of attack provided by the navigation 
subsystem and zeroing the filter rates 

^cmd _ OCnav' 

G cmd = a navi (initialization) (4) 

Vond = 0 . Am* = 0 

Commanded bank and angle of attack 
rates are generated, not by using the filter rates, but by 
comparing commanded angles with estimated angles and 
using a damping factor (1/2 here) to reduce overshoot. 

Anrf = («™, - a w )/ 2 

( 5 ) 

® cmd ~ (O cmd ~ ® nav) ^ ~ 

These rates are limited before being sent to the attitude 
control system 


The transition guidance algorithm is 
fairly simple. Desired entry interface attitude is 
specified via the mission design loads 
alphaTransitionDeg, and phibkTransitionDeg 


(6) 
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( 8 ) 



Heading Angles Definition 


u - {LI D)coser 

The horizontal distance (range) corresponding to flight 
along a drag-energy profile can be computed as follows. 
Define relative energy 

e = v 2 /2 + gh (9) 

and the time rate of change of energy is 


En try Guidance 

The guidance algorithm controls range to a 
specified landing site by issuing bank angle commands 
which will cause the reentry vehicle to track a nominal 
drag acceleration versus relative energy profile. The 
drag profile is a by-product of the overall (ascent 
through entry) trajectory design process, using the 
Program to Optimize Simulated Trajectories (POST). 
All relevant constraints such as range, heating and 
dynamic pressure limits are enforced in the trajectory 
design process. This ensures that when the reentry 
vehicle flies the resulting nominal drag profile, all 
relevant entry constraints are satisfied. Lateral trajectory 
control is provided via bank reversal logic which is 
described later. Entry guidance is terminated at an 
Earth-relative velocity of 2,500 ft/s. 

This section is organized as follows. 
First, range (longitudinal) control equations and logic 
are developed. Then, synthesis of a nonlinear tracking 
control law is discussed in which bank angle 
modulation is used to follow a reference drag-energy 
profile. Next, heading (lateral) control equations are 
discussed. 


e = vv + gh + gh - v(-D-gsiny) + gvsiny - - 

( 10 ) 


The time rate of change of range R is given by 
R = vcos y ~v (li) 


for near zero flight path angle (greater than about -5 
degrees). Therefore, the predicted range for flight on the 
profile can be computed from: 


dR „ 1 
de ~ D rtf 



( 12 ) 


In practice, if the drag profile is stored and used as a 
piecewise linear function of relative energy, (12) is 
integrable and the range calculation is a simple 
summation involving logarithms of the stored reference 
drag values. The summation can be performed a priori 
to obtain nominal range-to-go values corresponding to 
the tabulated energy values., reducing the nominal 
range-to-go calculation each guidance cycle to a simple 
table lookup. 


Range Control The basic idea is to develop a nominal 
drag versus relative energy profile which can be tracked 
by the vehicle via bank angle modulation. For guidance 
purposes, the relevant equations of motion are those 
corresponding to planar motion about a non rotating 
Earth, with central gravity and aerodynamic lift and drag 
forces 3 

r = vs'\ny (7) 

v = -D- gsin y 

V Y = ( v' 2 / r - g) cos y + Du 


Note that the restriction on flight path angle, i.e., small 
magnitude, implies that the valid region of reentry for 
drag-energy guidance extends down to about Mach 2.5, 
where typically flight path angle starts becoming large 
negative. Figure 3 shows a typical drag-energy profile. 

Tracking Control Law Synthesis In the following, the 
technique of feedback linearization is used to obtain an 
adaptive nonlinear tracking control law. Throughout 
the derivation of the control law, assume an exponential 
atmosphere model of the form 

p(h) = p 0 exp(-/3/i) (13) 


where the control u is the vertical or in-plane 
component of lift to drag ratio 


and use the approximation that time rate of change of 
drag coefficient is zero. (It is not difficult to remove the 
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latter assumption and doing so does not appreciably 
complicate the control law.) Define tracking error, Z » 
as the difference between (onboard) estimated drag 
acceleration and reference drag acceleration (from the 
stored drag profile): 

2(0 = D <A t )-D r , (0(14) 


z(t) = -/?£>„,( v si n(y) + vcos(y )y ) - (5rD„, + 

— v - 2-^v + 2^(-D t , - gsin(y) - 
v r’ v 

£cos(y)y - D' rt/ (D efl v + D. v ,v) 

= ~P(D es ,vs\n(y) + rD esl ) + —(£>„, - - 

V V 



-7 10 7 -6 10 7 -5 10 7 -4 10 7 -3 10 7 -2 10 7 -1 10 ? 0 

energy, fP/s 2 


Figure 3 

Typical Drag-Energy Profile 

A desirable time response for the tracking error is the 
solution of the following second order differential 
equation (the well-known damped harmonic oscillator) 
with specified values of damping ratio and natural 
frequency: 

z(t) + 2gw n z(t) + co-z(t) = 0 (15) 

In the following, expressions are obtained for the first 
and second time derivatives of tracking error. It will be 
seen that the time derivative of flight path (gamma dot) 
shows up (linearly) in the second derivative from which 
an expression for “commanded” gamma dot is obtained. 
From commanded gamma dot, the corresponding 
gamma dot dynamical equation is used to obtain a 
commanded in-plane lift component. First, z dot: 

on 

z(') = D,„ - D„ , = r + — «■ v - D' £>„,v 


The primes denote differentiation with 
respect to relative energy. The second time derivative 


-(D„,+gsin(y))- 


D' re/ (D esl v + D„,v) + Y(~j3D {5l vcos(Y ) - 2-^£cos(y)) 

V 


It is assumed that the drag-energy curve will be in the 
form of a linearly-interpolated table (i.e., D" e j = 0 ). 

Now the harmonic oscillator equation can be formed: 
-P(D esl vsm( Y ) + rD t J + — (£>„, ) - 

V V 

2D 

( D <s, + g sin(y )) - D' ref ( D es , v + D tu v) + y (-/?D„, v cos (y 
-2 — g cos(y)) + 2 £y(-)3D„ f r + ^££L v - D’ rtf D„, v) + 

V V 

DA = 0 


from which “commanded” gamma-dot is obtained: 


y cmd 

P( D tIl v sin(y ) + rD tsl ) + -(-£)„, + ^-) 

V V 

+— ( A, + ^in(y)) + D' rt/ (D esl v + D e J) - 

V 


2Cco(-pD es ,r + ^v 

V 


-D'D u v)-co;(D n ,-D rr/ ) 


~/?D„,vcos(y) - 2—^- gcos(Y) 

V 


Now this gamma-dot equation is equated with the 
gamma-dot equation (7) and solved for u, the inplane 
component of L/D: 

« = (vf , nu , - (V 2 / r - g)cos(Y))/ D e „ (20) 
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The control law depends upon the current slate (r, v, 
gamma), the drag-energy profile, the inverse scale 
height, beta, the estimated current drag acceleration and 
the design parameters omega and damping ratio, zeta. 
Note that the control law is independent of the current 
density and the current lift and drag coefficients. Thus, 
the control law has the interesting feature of not 
requiring a density model or an aerodynamic coefficient 
model. The estimated drag is just the component of 
sensed acceleration along the current navigated relative 
velocity vector. The sensed acceleration is equal to the 
accumulated velocity counts from the navigation system 
over the last one second. 


To compensate for downrange dispersions, two steps 
must be taken. First of all, the Terminal Area Energy 
Management (TAEM) target is effectively biased such 
that the current range error that exists is increased by 50 
percent. This is only used for control law evaluation 
purposes. This has the effect of driving the vehicle 
trajectory back to the nominal trajectory profile and 
preventing deviations from the nominal profile from 
growing as entry progresses. Secondly, the entire drag 
profile is (effectively) multiplied by the ratio of profile 
range-to-go to required range-to-go via: 


D _ R noJ e ) D r 


'ref 


R. 


r *f i nominal 


req 


D' - R non,( e ) p (<? ) 


ref 




nominal 


( 21 ) 


function of lift acceleration since it is the lift vector that 
is, in effect, being modulated to provide drag profile 
tracking. 

In addition to the basic control law above (a 
proportional-derivative control law), an integral control 
increment is computed 

A u=kA' (D„,-D ref )dt (23) 

W/ 0 

where the integration does not start until the transient 
dynamics have approximately died out, e.g. one period 
elapsed since closed-loop guidance initiation. The 
intent of the integral portion of the control law is to 
eliminate steady state error. 5,67,8,9 The gain K on the 
integral term should be chosen very small to preserve 
the desired transient response properties as specified 
with the damping ratio and natural frequency. For the 
nonlinear system under consideration, engineering 
judgment and high-fidelity simulation testing are 
required to ensure that the integral portion of the control 
law enhances tracking performance and stability rather 
than degrades it. 

Angle of attack is nominally commanded according to a 
pre-specifted Mach-alpha schedule. Figure 4 shows a 
typical angle of attack profile. To augment bank angle 
control, angle of attack modulation about the nominal 
schedule is used. Assuming drag coefficient is linear 
with alpha, an incremental alpha, A CC r required to 
make current drag equal to desired drag is 


where R nom is the nominal range-to-go computed using 

the nominal drag profile. In the presence of heading 
error, a reasonable approximation of required range-to-go 
is given by (for a circle on a flat earth 4 ): 


R. 


req 


s% 

sinOF,) 


( 22 ) 


D ,-D , = A D = 

ref est 

pv 2 5AC D ^ pv 2 5C Po Aa 
2m 2m 


=> A a = 


ADC d 


P M C Dm Aa (24) 

Q> 

(25) 


where x ¥ e is the current heading error The cyclically 

calculated values from (21) are the values of reference 
drag and slope that are used in the control law. 

A value of 0.7 is highly recommended for the damping 
ratio 4,5,6,7,8 In Ref. 6 the period (inversely proportional 
to CO n ) is a function of drag acceleration level. A high 

value of drag implies that more aerodynamic force is 
available, which means that a faster response (smaller 
period) is possible for a given amount of bank 
maneuvering. In X-33 guidance, the period is a linear 


This value is passed through a limiter to keep 
commanded angle of attack within five degrees of the 
nominal schedule. This alpha modulation capability is 
particularly useful during roll reversals and during short- 
period transients. An adverse effect of alpha 
modulation is that it's effect can sometimes be the 
opposite of what is desired. For example, to decrease 
drag, a lower value of alpha is commanded which also 
reduces lift which reduces altitude which increases drag. 
This adverse effect is especially noticeable when bank 
angle is saturated at all lift up. 
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A bank angle increment, the intent of 
which is to drive the angle of attack profile back to the 
nominal schedule, is added to bank command 6 : 

= cos' 1 ((« + A u)Dt L) + 4A a (26) 

The bank command is limited to between 0 and 85 
degrees so that when a roll reversal is commanded, 
smaller maneuvers are required. The sign of the bank 
command is determined from the logic of the next 
section. Finally, to ensure that bank and angle of attack 
commands are reasonably consistent with vehicle 
attitude maneuver limitations, these two quantities are 
passed through a function that models the bank and 
angle of attack dynamics as second order systems with 
pre-specified rate and acceleration limits. This is 
analogous to the Transition Guidance logic. As in 
transition guidance, outputs to the control system are 
again bank, bank rate, alpha, and alpha rate. 

Heading Angle Control The lateral logic is, for the 
most part, uncoupled from the longitudinal (range) 
control. During reentry, the lateral logic determines the 
sign of the commanded bank angle and does it 
independently of the range control (which determines the 
appropriate magnitude of the bank angle command). 

In Ref. 6, a heading error versus speed corridor is stored 
onboard to determine when to command roll reversals. 
To minimize the extent to which the heading error 
corridor is tied to a specific entry profile, a new, more 
general approach has been developed for X33 reentry 
lateral control, a description of which follows. 

The time rate of change of velocity heading angle is 
written 3 

Lsin(cr) v 

¥ = — + - cos(/)sin(yf )tan(0) ~ 

vcos(y) r (27) 

2to(tan(x)cos(y/)cos(0) - sin(0» 

Likewise, an expression for the time rate of change of 
desired heading (omitted here for brevity) can be 
obtained from its definition so that the time rate of 
change of heading error is 


Lsin(cr) 

= V (i -7 + 

vcos(x) 

V' 

-cos(x)sin(i//)tan(0) - (28) 

r 

2w(tan(/)cos(^)cos(0) - sin (<p)) 

Equating to zero, and solving for bank angle gives 
2 

sin(a) = \j/ d + — cos 2 (y)sin(>)tan(0)- 

Lr 

, . (29) 

_ vcos(y) 

2(0 —= — (tan(y)cos(i//)cos(0)-sin(0)) 

Lj 

This value is the bank angle required to maintain the 
current heading error. If the magnitude of the required 
bank is greater than 12 degrees, then a sign change (on 
the commanded bank angle) is indicated. 

When Mach decreases to below 4.0, bank angle is no 
longer determined from range control considerations but 
rather from heading consideration, i.e., bank angle 
command is proportional to the current heading error. 

a = 4 (v d - v) (30) 


Performance Monitor 

On-board the X-33, a separate computer contains the 
Mission Manager (MM) logic. The Mission manager 
detects the current phase of flight and executes the on- 
board Performance Monitor (PM) which is a 3 degree- 
of-ffeedom X-33 simulation that predicts the trajectory 
from the current vehicle state to the TAEM interface. 
Under certain conditions, the PM may send commands 
to the FM which modifies the nominal guidance 
steering command to the control system. In these 
situations the MM/PM functions as a closed-loop 
guidance algorithm. 

One of the duties of the Mission Manager is to detect 
the depletion of propellant and signal a change of flight 
mode fron the ascent to transition phase of flight. 
Another important role of the Mission Manager is to 
detect an engine out situation and, based on the time of 
failure, use the Performance Monitor to reshape the 
trajectory to the nominal landing site, reset the guidance 
targets to divert to an alternate landing, or predict the 
need for an intentional destruction of the vehicle. 
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The Performance Monitor is a 3- DOF simulation that is 
intended to be flown on-board the X-33 vehicle as pan 
of (light software. The Performance Monitor will be 
executed continuously, initializing to current flight state 
data and executing until the TAEM target conditions are 
either missed or achieved. These results will provide a 
prediction, in real time, whether a safe landing at the 
nominal landing site can be achieved. The Performance 
Monitor will also be used to assess the reshaped 
trajectory for an open-loop ascent abort before allowing 
the guidance system to use the new real-time computed 
chi tables. The on-board trajectory reshaping problem is 
complicated by a minimum acceleration constraint on 
the engine to avoid pump cavitation. The minimum 
acceptable axial acceleration is a function of the 
remaining liquid oxygen and a function of the throttle 
level. 

The MM uses a drag factor parameter as an indicator of 
how much energy margin the vehicle will have at the 
beginning of the entry phase. The drag factor is 
calculated within the PM and is the ratio of the drag 
profile range to the required range to reach the nominal 
TAEM conditions. The drag profile range is the 
downrange distance that the vehicle could achieve if at 
the current energy condition , the vehicle followed the 
nominal drag profile. The required range is the great 
circle range from the current vehicle position to the 
TAEM point. The drag factor decreases significantly for 
early engine pump failure times because the energy 
level at MECO becomes much lower. 

Due to the decreased vehicle thrust to weight ratio if one 
of the two linear aerospike engines fails early in flight 
it is not always possible to reach the nominal or 
alternate landing site. In this case the Mission Manager 
would predict the need to destroy the vehicle. 
Fortunately, the vehicle spends a considerable amount 
of time over the EAFB site so the vehicle impact point 
would be well within EAFB boundaries. The minimum 
throttle range constraint does not allow’ for RTLS-type 
trajectories. 

Figure 5 shows the final range to HAC as a function of 
drag factor. It was found that when the drag factor is 
above 0.7, the closed loop entry guidance could 
successfully achieve TAEM conditions without any 
reshaping of the ascent trajectory by the Performance 
Monitor. This corresponds to single engine failure 
times greater than 40 seconds and all of the likely 
vehicle and trajectory dispersions that have been studied. 

The trajectory profiles in Figure 6 indicate that the 
trajectories for failure times greater than 40 seconds arc 


lofted. This lofting is induced because at a given 
altitude, the vehicle has less velocity than it would on a 
nominal trajectory. Since the ascent open loop attitude 
commands (Chi Table) are based on velocity, the pitch 
commands at a given altitude are higher than in the 
nominal case. For late engine failure times entry 
guidance can compensate for the lower energy state at 
MECO. However, for early failure there is not 
suffucient energy to achieve the TAEM targets unless 
the ascent trajectory is modified. 

The trajectory profiles for early engine failures show the 
effect of Performance Monitor reshaping on the ascent 
trajectories. Initially, just after the pump failure, the 
Performance Monitor reshaping logic reduced the value 
of the pitch command from the open loop Chi table. 
This reduces the trajetory lofting. During the later 
portion of the ascent trajectory the value of the pitch 
command is increased. This modification (or reshaping) 
of the pitch commands in the nominal Chi table enables 
the vehicle to reach an acceptable energy level at MECO 
which allows entry guidance to achieve the target 
TAEM conditions. 

The use of an on-board 3-DOF vehicle simulation to 
evaluate, in real time, the ability of the guidance 
system to achieve its intended target objective has not 
been attempted before on any launch vehicle or 
spacecraft. This innovative X-33 flight software design 
will add significant “intelligence” to the GN&C system. 

Dispersion Results 

The ascent, transition and entry guidance algorithms 
were tested in a 3 DOF simulation against 
environmental and system dispersions. The various 
parameters were dispersed in a Monte Carlo fashion. 
Dispersion studies were ran for three of the five test 
flights to MAAF. The dispersion results to be 
presented here are from the third flight: Michael 7c3. 

This purpose of this particular flight is to reach a 
specific level of catalytic aeroheating. Table 1 lists the 
parameters which were dispersed and the level of their 
dispersions. An EAFB 3hr wind-pair database was used 
for this dispersion analysis. A wind-biased trajectory 
was designed for the first of the wind and simulated 
against the second. The GRAM 97 (Global Reference 
Atmosphere Model) was employed to randomly vary the 
atmosphere for a given month and the winds during 
entry. Results here will be limited to the set of January 
simulations. 
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Dispersion 

Parameters 

Dispersion 

Level 

Aerospike Engine 


Mixture Ratio Error 
(LOX:LH2) 

± 2% 


± 1% 

LOX Flow Rate Error 

± 0.6% 

LH2 Flow Rate Error 

± 1.125% 

Environments 


Ascent Atmosphere 

GRAM 97 EAFB 
Jan, July & Nov 

Entry Atmosphere 

GRAM 97 
Jan, July & Nov 

Ascent Winds 

2nd of 
EAFB DOL 

Entry Winds 

Randomly 
Dispersed Wind 
Profile 

Mass Properties 



± 0.6% 

Loaded LOX Error 

± 0.232% 

Loaded LH2 Error 

± 0.181% 



CA/CN 

5% 

Navigation 

Errors 

Navigation 
uncertainty model 
specified in 
DFRC Dispersion 
Document vi.l 


Iflfels 1 

Monte Carlo Dispersion Parameters 

Many output parameters were tracked during each of the 
Monte Carlo simulations. Ascent loading indicators, 
MECO parameters such as cutoff time, altitude, Mach 
number, and flight path angle. Entry parameters were 
also tracked including thermal and structural indicators. 
The results presented in this paper will be limited to the 
parameters recorded at the TAEM interface (defined when 
the relative velocity reaches 2500 ft/s). In this manner, 
we can judge how well the ascent, transition and entry 
guidance steered the vehicle to the final target.. The 
tracked TAEM parameters are listed in table 2. 

The guidance performed well under the given dispersion 
levels. There were extreme cases that required further 
investigation but the standard deviation values were 
with the TAEM target zone. 

Figure 7 is a plot of the maximum and minimum 
crossrange as a function of the vehicle downrange for 


this set of disperison simulations. Again, the 
groundtrack spread was within acceptable values. 


Variable 
Taem Time 

altFT 

Mach 

fpaDeg 

rangeLand 


Description 

Trajectory time at which the 
TAEM condition is reached (s) 
Altitude at TAEM (ft) 

Mach number at TAEM 
Flight path angle at TAEM (deg) 
range to the TAEM point (nmi) 


Tzbk Z 

Track ed TAEM Di s pe rs io n Pa r ame ters 



Min 

Max 


Stdev 

Taem 

Time 

428 

464 

445.6 

5.5 

altFt 

84394.1 

101439 

94018.4 

2044.6 

Mach 

2.44 

2.68 


0.05 


0.1 

19.9 

1.67 

1.44 


-17.8 

-6.9 

-11.2 

1.4 


Table 3 

Michael 7c3 Trajectory. Dispersion Results 
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Figure 4 

Typical Entry Angle of Attack Profile 


Range to HAC as a Function of Drag Factor for the 
Michael 7c7 X-33 Mission 



Figure 5 

Range to HAC as a Function of Drag Factor 
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X-33 Power Pack Failure Guided Abort Trajectories tc 
Michael 

(With Performance Monitor Reshaping) 



Rang* to Landing Site (nml 


Figure 6 

X-33 Nominal and Abort Trajectory Profiles 


Michael 7c3 January 
(No PM reshaping or PU system model) 



Downrange (nm) 

Figure 7 

X-33 Trajectory Dispersion Results 
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